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Abstract. We study the warping and tearing of a geometrically thin, non-self-gravitating 
disk surrounding binary supermassive black holes on an eccentric orbit. The circumbinary disk 
is significantly misaligned with the binary orbital plane, and is subject to the time-dependent 
tidal torques. In principle, such a disk is warped and precesses, and is torn into mutually 
misaligned rings in the region, where the tidal precession torques are stronger than the local 
viscous torques. We derive the tidal-warp and tearing radii of the misaligned circumbinary 
disks around eccentric SMBH binaries. We find that in disks with the viscosity parameter 
a larger than a critical value depending on the disk aspect ratio, the disk warping appears 
outside the tearing radius. This condition is expressed as a > for H/r < 0.1, 

where H is the disk scale height. If a < only the disk tearing occurs because the 

tidal warp radius is inside the tearing radius, where most of disk material is likely to rapidly 
accrete onto SMBHs. In warped and torn disks, both the tidal-warp and the tearing radii 
most strongly depend on the binary semi-major axis, although they also mildly depend on 
the other orbital and disk parameters. This strong dependence enables us to estimate the 
semi-major axis, once the tidal warp or tearing radius is determined observationally: For 
the tidal warp radius of 0.1 pc, the semi-major axis is estimated to be 10 ^ pc for 10^ Mq 
black hole with typical orbital and disk parameters. We also briefly discuss the possibility 
that central objects of observed warped maser disks in active galactic nuclei are supermassive 
black hole binaries. 


^Corresponding author. 
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1 Introduction 

Supermassive black holes (SMBHs) with mass 10® Mq < M < are considered to 

reside at the center of most galaxies [1]. Hitherto, SMBHs have been found in 87 galaxies by 
observing the proper motion of stars bound by the SMBHs or by detecting radiation emitted 
from gas pulled gravitationally by the SMBHs [2]. H 2 O maser emission from active galactic 
nuclei (AGNs) in spiral galaxies provides a strong tool to measure supermassive black hole 
(SMBH) masses, because it arises from a rotating disk on a subparsec scale with a nearly 
Keplerian velocity distribution around a SMBH. Those maser disks have been observed at the 
centers of NGG4258 [3], NGG 1068 [4], NGG3079 [5], the Gircinus galaxy [6], UGG3789 [7], 
NGG6323 [8], NGG2273, NGG6264, and some more objects [9]. Several maser disks show 
warped structure at the radii of ~ 0.1 pc [6, 9-11]. From an observational point of view, 
maser spots on the disk in NGG 4258 are spatially distributed along a line on each side of a 
central black hole. The SMBH is then thought to be located at the center of a line connecting 
those two lines by a simple extrapolation, and the disk starts to warp at the innermost maser 
spot. What mechanism makes the disk warped still remains an open question. 

Several scenarios have been proposed for explaining disk warping. Pringle [12] showed 
that centrally illuminated accretion disks are unstable to warping due to the reaction force 
of reradiated radiation. Such a radiation-driven warping mechanism has also been applied to 
explain the disk warping in the context of X-ray binaries [13-16]. If angular momentum vector 
of an accretion disk around a spinning black hole is misaligned with the spin axis, differential 
Lense-Thirring torque due to the frame-dragging effect aligns the inner part of the disk with 
the black-hole equatorial plane. Since the outer part of the disk retains its initial orientation, 
the resultant disk is warped [17]. This Bardeen-Petterson effect is also considered to be a 
plausible mechanism for disk warping in maser disks [18]. Moreover, Bregman & Alexander 
[19] proposed that the warped disk at the center of NGG 4258 is caused by the process of 
resonant relaxation, which is a rapid relaxation mechanism to exchange angular momentum 
between the disk and the stars moving under the nearly spherical potential dominated by 
the SMBH. These mechanisms have been mostly discussed based on the assumption that the 
central object surrounded by the warped maser disk is a single SMBH. 

The tight correlation between the mass of SMBHs and the mass or luminosity of the 
bulge of their host galaxies ([20-23]; see also [2] for a review) suggests that the SMBH at the 
center of each galaxy should have evolved toward coalescence in a merged galaxy. If this is 
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the case, a binary of SMBHs on a parsec or subparsec scale should be formed in a merged 
galactic nucleus before two black holes finally coalesce, yet no SMBH binaries have clearly 
been identified so far despite some claims (see [24, 25] for reviews and references therein). 

A possible link between the presence of SMBH binaries and the warping of observed 
maser disks has been studied by Hayasaki et al. [26] (hereafter, H14a). They have, for 
simplicity, assumed that the circumbinary disk is initially aligned with the binary orbital 
plane. However, the angular momentum vector of the circumbinary disk does not always 
coincide with that of the binary orbital angular momentum, because the orientation of the 
circumbinary disk is primarily due to the angular momentum distribution of the gas supplied 
to the central region of AGNs. Therefore, the orientation of the circumbinary disk plane is 
expected to be distributed randomly with respect to the binary orbital plane. The formation of 
such a misaligned circumbinary disk around SMBH binaries has been numerically examined by 
Dunhill et al. [27]. Hayasaki et al. [28] have investigated the accretion process from misaligned 
circumbinary disks onto SMBHs in eccentric orbits. 

Regarding the misaligned disk structure, the inner part of the circumbinary disk tends 
to align with the binary orbital plane, because the tidal precession torque is stronger than 
the shear viscous torque in the vertical direction, whereas the outer part tends to retain the 
original state because the tidal precession torque is weaker than the vertical viscous torque. 
This is the origin of disk warping [29-31]. 

On the other hand, Nixon et al. [32] have recently proposed that the circumbinary disk 
is broken into mutually misaligned rings, if the tidal precession torque is stronger than the 
local horizontal viscous torque in circular SMBH binaries. However, little is known about the 
relationship between disk tearing and warping in eccentric SMBH binaries. 

In this paper, we examine the tidally driven warping and tearing of misaligned cir¬ 
cumbinary disks around eccentric SMBH binaries. In Section 2, we describe the tidal torques 
originating from a time-dependent binary potential and derive the tidal warp radius of the 
misaligned circumbinary disk. In Section 3, we discuss a possibility that the observed warp¬ 
ing of maser disks in several AGNs is caused by the tidal effect of SMBH binaries. Finally, 
Section 4 summarises our scenario. 


2 Tidally driven warping and tearing of a misaligned circumbinary disk 

Let us consider the torques from the binary potential acting on the circumbinary disk, which 
is misaligned with the binary orbital plane, surrounding two black holes in a binary on a 
eccentric orbit. Figure 1 illustrates a schematic picture of the setting of our model; binary 
black holes orbiting each other are surrounded by a misaligned circumbinary disk. The 
binary is put on the x-y plane with its center of mass being at the origin in the Cartesian 
coordinate. The masses of the primary and secondary black holes are represented by Mi and 
M 2 , respectively, and M = Mi -|- M 2 . We put a circumbinary disk around the origin. The 
unit vector of specific angular momentum of the disk is expressed by (e.g. Pringle [12]) 

I = cos 7 sin /3x -|- sin 7 sin j3y + cos (3z, ( 2 . 1 ) 

where /3 is the tilt angle between the circumbinary disk plane and the binary orbital plane, 
and 7 is the azimuth of tilt. Here, x, y, and z are unit vectors in the x, y, and z, respectively. 
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Table 1. Model parameters 


Dehnition 

Symbol 

Total black hole mass 

M 

Primary black hole mass 

Ml 

Secondary black hole mass 

M 2 

Schwarzschild radius 

rs = 

Binary mass ratio 

q = M 2 /Ml 

Mass ratio parameters 

6 = 9/(1 + 9), 6 = 1/(1 + g) 

Binary semi-major axis 

a 

Orbital eccentricity 

e 

Orbital frequency 

Oorb = y/GM/a^ 

Orbital period 

Porh — 27 r/norb 

True anomaly 

+ 

II 

Tilt angle 

/3 

Azimuth of tilt 

7 

Azimuthal angle 


Shakura-Sunyaev viscosity parameter 

a 

Horizontal shear viscosity 


Vertical shear viscosity 


Ratio of vertical to horizontal shear viscosities 

V = ^21^1 

Mass-to-energy conversion efficiency 

e 


The position vector of the disk can be expressed by 

r = r(cos(/>sin 7 + sin(/)cos 7 cos/ 3 )a: 

+ r (sin (j) sin 7 cos /? — cos 4> cos 'y)y 

— r sin (/) sin( 2 . 2 ) 


where the azimuthal angle (p is measured from the descending node. The difference from 
equation (3) of Paper I is the position vector of each black hole, which is given by 

ri = riCosfiX + risii\fiy (z = l,2), (2.3) 


where /2 


/i + TT is the true anomaly and r* is written as 


ri = 


0(1 — e^) 

1 + e cos fi 


(2.4) 


with ,^1 = g/(l + ( 7 ) and ^2 = 1/(1 +?)• Here, e is the binary orbital eccentricity, q = M 2 /M 1 
is the binary mass ratio, and a = ai + 02 is the binary semi-major axis with ai = and 
02 = ^ 20 - These and other model parameters are listed in Table 1. We also assume that the 
disk total angular momentum is small compared with the binary angular momentum. 


2.1 Tidal torques acting on the misaligned circumbinary disk 

The gravitational force on the unit mass at position r on the circumbinary disk can be written 

by 


F 


grav — 


E 


GMi 



{r 


ri) 


(2.5) 


3 








Tilt angle:/3 



Figure 1. Configuration of a SMBH binary on an eccentric orbit and a circumbinary disk surrounding 
it. There are two angles (/3, 7 ) which specify the orientation of the disk plane with respect to the 
binary orbital plane (x-y plane). The azimuthal angle (</)) of an arbitrary position on the disk is 
measured from the descending node. 


The corresponding torque is given by 


t 


grav 


2 

V X F grav — ^ ^ 
i=l 


GM, 



(r X Vi) 


( 2 . 6 ) 


We consider the tidal warping/precession with timescales much longer than local rotation 
period of the circumbinary disk. This allows us to use the torque averaged in the azimuthal 
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direction and over the orbital period: 


p27T i'2n 

(Tgrav) ^ ^^2 J J ^grav 

_ 1 r^^ll ^GMi{r xvi) ^2 
2vr Jo 27r io ^ 


GMi{r X rj) 
T3— 


r - Vi 


(1 + ecos/i )2 


df 


3 . . GM /a \2 

= -( - 

8 r \r 


+ 5e^ sin 2'y sin^ f3z 


(1 — e^) sin 7 sin 2/3® — (1 + 4e^) cos 7 sin 2/3y 


(2.7) 


which is equivalent to equation (7) of Hayasaki et al. [33], where florb = \/GM /is the 
angular frequency of the mean binary motion. Here, we used for integration the following 
relationship: 


(I _ e 2 ') 3/2 

^ (l + ecos/O^"^-^ 


and the approximations: 


r — r, 


1-3 


„-3 


l + 3^ + 0 ((r,/r)' 


(1 + ecos fi) « 1 — 4ecos fi + lOe^ cos"^ fi + 0{e'^), 


( 2 . 8 ) 


(2.9) 

( 2 . 10 ) 


where r S> a is adopted for equations (2.9). The magnitude of the specific tidal torque is 
given by 

Krg„.)l = (^)"e(e./S.7) (2.11) 

where 0 (e,/ 3 , 7 ) is a function given by 


0(e,/3,7) = { 4 


5e^(2 + 3e^) cos^ 7 + (e^ — 1)^ 


cos^ /3 


1/2 


+ 25e^ sin^ 27 sin^ /3 ^ sin /3. 


( 2 . 12 ) 


Figure 2 shows dependence of 0(e,/ 3 , 7 ) on the three parameters: orbital eccentricity e, tilt 
angle /3, and azimuth of tilt 7 . It is noted from the figure that there are several zeros in 
0(e, /3, 7 ), e.g., for /3 = 0 and vr with any e and 7 . Except for these values, 0(e, /3, 7 ) is of the 
order of one. 

By going through the same procedure as for equation (2.7), we obtain the azimuthally- 
averaged and orbit-averaged tidal potential acting on the misaligned circumbinary disk around 
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Figure 2. Dependence of 0(e, /3, 7) on orbital eccentricity e, tilt angle /?, and azimnth of tilt 7. Panel 
(a); dependence on e with fixed valnes of /3 and 7. The solid black, red dashed, and bine dotted lines 
indicate 0(e, 7 r/ 4 , 7 r/ 4 ), 0(e, 7 r/ 2 , 7 r/ 4 ), and 0(e, 7 r/ 4 , 7 r/ 2 ), respectively. Panel (b); dependence on (3 
with fixed valnes of e and 7. The solid black, red dashed, and bine dotted lines indicate 0(0.0, (3, 7 r/ 4 ), 
0(0.6,/3, 7 r/ 4 ), and 0(0.9,/3, 7 r/ 4 ), respectively. Panel (c); dependence on 7 with fixed valnes of e and 
(3. The solid black, red dashed, and bine dotted lines indicate 0(0.0, 7 r/ 4 , 7 ), 0(0.6, 7 r/ 4 , 7 ), and 
0 ( 0 . 9 , 71 / 4 , 7 ), respectively. 


an eccentric binary: 


$ = 


1 A GMi 




rZ7T /•ZTT ^ 

LI 


GM 


1 + 


66 


2=1 

. 3 2 


\r - Ti 


-d(f)d{yt,orht) 


a\ 2 
r ) 


366 


(1 — + 5e^ cos^ 7 ) sin^ (3 


(2.13) 


Here, we used for integration equation (2.10) and the following approximation: 


\r — rA ^ Ri r ^ 


1 + 5^7)+ ^( 35 ^- 1 ) (7) +^(666^) 


(2.14) 


where gi = cos (/)sin (7 — /j) + sin (/> 003(7 ~ fi) cos (3. Since the tidal potential has minima at 
/3 = 0 (alignment) and (3 = n (counter-alignment) for arbitrary e and 7 , the tidal torques 
tend to align or counter-align the tilted circumbinary disk with the orbital plane. Note that 
equation (2.13) is reduced to equation (3) of Hayasaki & Okazaki [34] if ,0 = 0. 

The tidal precession timescale for eccentric binaries is given by 

1 

66 

_ 4 / 1 \ /r\7/2 1 

3vr V66y 0 (e, 7 ,/ 3 ) 

where J = r^Ql with the disk angular frequency D = GM/r^ and Porb = 27r/Dorb are 
the specific angular momentum and binary orbital period, respectively. Note that the tidal 
precession timescale depends on the tilt angle and azimuth of tilt for e 7 0- Since the inner 
edge of the disk is estimated to be ~ 2a [35] , the tidal precession timescale is longer than the 
binary orbital period. 


r \ 2 


a/ 0(e, 7 , /3) n 


(2.15) 
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2.2 Tidal warping vs. tearing 


There are two types of important viscosities in the circumbinary disk. The first type is the 
horizontal shear viscosity, ui, which is the viscosity normally associated with accretion disks. 
The second type is the vertical shear viscosity, V 2 , which tends to smooth out disk warping 
when the disk is non-planar. Ogilvie [36] derived the relationship between vi and V 2 '. 


^ ^ 2(1 + 7a^) 

vi 


(2.16) 


by taking a non-linear effect of the fluid on the warped disk, where r/ and a are the viscosity 
ratio parameter and the Shakura-Sunyaev viscosity parameter, respectively. For a <C 1, the 
above equation is reduced to 


7] ^ 


1 


(2.17) 


The global horizontal and vertical viscous timescales for a geometrically thin disk are 
estimated to be 




I 1/2 


2r2 ^ 

^211/ 

3 ui 

" SaQ V 

2r2 ^ 

^211 

3 1^2 





(2.18) 

(2.19) 


where i^i ~ aClH^ with the disk scale hight H and 1/2 ~ rjaQH^. On the other hand, the 
local horizontal viscous timescale is estimated to be 


_ 2rAr _ 2 1 1 fH\~^ 

3 ui 3all\ry ’ 


( 2 . 20 ) 


where Ar ~ 77 is adopted. 

The misaligned circumbinary disk is warped at a radius where the tidal precession 
timescale is longer than the vertical viscous timescale. We call this radius the tidal warp 
radius, following Martin et al. [14, 15], who used it in the context of X-ray binaries. The tidal 
warp radius is given by 


'warp 

rg 


3 66 

8 r]a 


1/2 


|0(e,/3,7)| 




( 2 . 21 ) 


Recently, Nixon et al. [32] have proposed that the circumbinary disk is torn into mutually 
misaligned gas rings if the timescale of the tidal precession, which is strongly differential, is 
shorter than the local horizontal viscous timescale. The disk tearing occurs at the radius: 


72tear 

rs 


1 66 
4 a 


1/2 


|0(e,/3,7)| 




( 2 . 22 ) 


This reduces to equation (9) of Nixon et al. [32] when e = 0. The disk tearing is also 
conhrmed by a three-dimensional Smoothed Particle Hydrodynamics simulations [32]. Inside 
the tearing radius, the material rapidly accretes onto the central binary, especially, in a 
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Figure 3. Critical value of viscosity parameter as a function of H/r. The black solid line is the 
critical value, Uc, of a without approximation to 77, whereas the red dashed line is ac for rj Ri 1/(2 q!^). 

retrograde, misaligned circumbinary disk. Therefore, in order for the disk to have a warped 
structure, i?warp must be larger than iJtear- This condition is equivalent to 

3 1 H 
2rj r 

This is then rewritten as 

a > ac, (2.24) 

where ac represents the critical value for the disk viscosity parameter, which is a function of 
H/r. Figure 3 shows the dependence of ac on H/r. The red dashed line denotes the critical 
value of a for a <C 1, which is equal to y/H/{3r) obtained by substituting equation (2.17) 
into equation (2.23). From the figure, we note ac = \/H/{3r) is a good approximation for 
geometrically thin disks with H/r < 0.1. For a > etc, the disk warping occurs outside the 
tearing radius. Otherwise, no disk warping will occur because there is little disk material at 
the tidal warp radius, which is inside the tearing radius. 

Figure 4 shows the dependence of tidal-warp and tearing radii on the semi-major axis. 
For a fixed value of H/r = 0.01, ac ~ 0.06. We can conhrm that the tidal warp radius is 
larger than the tearing radius for any semi-major axis for a = 0.1, whereas it is smaller than 
the tearing radius for any semi-major axis for a = 0.01. 

The disk evolves diffusively for a > H/r, whereas it does with dispersive wave propaga¬ 
tion for a < H/r [37]. For the latter case, the communication time due to wave propagation 
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Figure 4. Dependence of tidal-warp and tearing radii on the semi-major axis. The black solid 
and dashed lines are the tidal-warp and tearing radii for a > ac case, respectively, whereas the 
red solid and dashed lines are the tidal-warp and tearing radii for a < Uc case, respectively. Here, 
H/r = 0.01 (oc 0.06) is adopted for each radins. The bine line shows the inner edge radins of the 
circnmbinary disk. The tidal-warp and tearing radii have other common parameters: M = IO^Mq, 
q = 0.1 (^ 1^2 = 10/121), e = 0.6, j3 = 7 r/ 4 , and 7 = 7 r/ 4 , respectively. 

is roughly given by 2r where Cg is the local sound speed. The tearing radius is then given 

by 


R. 


tear,w 

r-g 


^66|0(e,/3,7)l 


1/2 


r ) I rs 


Comparing this equation with equation (2.21), we obtain, for a <C 1, 


R 


warp 


( 4 ) 


r \i/2 


< 1 . 


-Rtear,w 

Thus, no disk warping occurs in a tilted circmbinary disk with a < H/r. 


(2.25) 


(2.26) 


2.3 Warping and tearing of circnmbinary disks in AGNs 

In this section, we apply our model to AGN disks. A gaseous disk around a SMBH in an 
AGN is surrounded by a dusty torus (e.g. [38]). The grains in the dusty torus is evapo¬ 
rated above the dust sublimation temperature by the irradiation emitted from the central 
source. The inner radius of the dusty torus should therefore be determined by the dust 
sublimation radius: r^ust = 3pc (L/10^® ergs“^)^/^(Tdust/lh00K)“^'®, where T^ust is the 
dust sublimation temperature [39]. Assuming that the AGN luminosity is the Eddington 
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luminosity: LEdd — 1-3 x 10 ^®(M/Mq) ergs“^, the dust sublimation radius is rewritten as 
^dust = 4.7 X 10“^ (M/10^ pc with the adoption of Tdust = 1500 K. Since the cir- 

cumbinary disk should be also inside the dusty torus in our scenario, we assume the disk 
outer radius and temperature at iiout to be given by 

«4.8 X 10^ ( ' (2.27) 

and Tout = 1500 K, respectively. 

Here, we assume that the radial profile of the disk temperature obeys 


T — Tout 



(s < 0), 


(2.28) 


where s, Tout) and Tout are the power law index of the radial profile of the disk temperature, 
disk outer radius, and temperature at Tout: respectively. With the thin disk approximation, 
the ratio of the disk scale hight to radius is then given by 


H 

r 



(2.29) 


where Cg^out = \/ (Tg//u)T is the isothermal sound speed at Tout with the gas constant Tg 
and molecular weight /r. 

Substituting equation (2.29) into equation (2.21), the tidal warp radius is rewritten as 


T 


'warp,AGN 

r-g 


3 6^2^, O 

- 0 (e,/ 3 , 7 ) 


16 rja 


CsjOut 


2/(s+3) 


T, 


out 


rg 


s/(s+3) 


-) 


2/(s+3) 


(2.30) 


Similarly, the tearing radius is rewritten as 
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tear,AGN 

rg 


1 6‘^2 


_4\/2 a 


0 (e,/ 3 , 7 ) 


T, 


out 


rg 


s/(s+5) 


-) 


2/(s+5) 


4/(s+5) 


Cs.out 


2/(s+5) 


(2.31) 


Figure 5 shows characteristic radii of the misaligned circumbinary disk for e = 0.6, 
q = 0.1 (Ci ?2 = 10/121), M = IO^Mq, a = 0.1 {rj k, 50), s = —3/4, /3 = 7 r/ 4 , and 7 = 7 r/ 4 . 
The solid and dashed black lines show the tidal warp and tearing radii, respectively. The red 
solid line shows the radius, which we call the orbital decay radius, where the tidal precession 
timescale equals the timescale in which the binary orbit decays by the gravitational wave 
(GW) emission. This orbital decay timescale for an eccentric binary is given by [40] 


5 1 


Tgw — 


’ 66 


rs (1 - e2)7/2 


c l + 73e2/24 + 37e4/96 
Equating equation (2.32) with equation (2.15), we obtain the orbital decay radius: 

R , r 7/2 1 \ 12/7 


prec/gw 


15 0 (e,/ 3 , 7 )(l - 

6462 1 + 73eV24 + 376^/96 


-) 

rsj 


(2.32) 


(2.33) 
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Figure 5. Characteristic radii of the warped and torn circumbinary disk around the SMBH binary 
on an eccentric orbit for e = 0.6, a = 0.1 {r] k. 50), s = —3/4, q = 0.1 (^ 1^2 = 10/121), M = IO^Mq, 
(3 = 7 r/ 4 , and 7 = 7 r/ 4 . The solid and dashed black lines show the tidal warp and tearing radii, 
respectively. The red solid line shows the orbital decay radius where the tidal precession timescale 
equals the orbital decay timescale due to the GW emission (see equation (2.33)). While the blue lower 
line represents the inner radius of the circumbinary disk rin ~ 2 a, the blue upper line represents the 
outer radius of the circumbinary disk i?out/?'s ~ 4.8 x 10"^ (M/10^ . The shaded area between 

the two blue lines represents the whole region of the circumbinary disk. 


When it is larger than the tidal warp and tearing radii, the circumbinary disk can be warped 
and torn before two SMBHs coalesce. Otherwise, two SMBHs coalesce prior to disk warping 
and tearing. The blue solid lower and upper lines show the inner and outer radii of the 
circumbinary disk, respectively. The shaded area between the two blue lines shows the whole 
region of the circumbinary disk. 


3 Discussion 


In this section, we discuss a possible link between the observational signatures of circumbinary 
disk warping and tearing and the presence of SMBH binaries. If the warp and/or tearing radii 
are observationally determined, we can estimate the semi-major axis of the SMBH binary from 
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equations (2.30) and (2.31) as 
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®tear 
rg 
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16 f r]a 
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1 1/2 


3 V 66 y 0(e,/3,7) 
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Os,out \ / .otout \ 
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/ ^obs \ 

V I’s y 


(s+3)/2 


4^2 
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1 1/2 


/ ^obs \ 

V I’s y 


66 y 0(e,/3,7) 

(s+3)/4 


Cs,out\V2 


rg 


(3.1) 


(3.2) 


where i?obs shows the tidal warp or tearing radius determined by observations. Figure 6 
shows the relationship between this observed radius of disk warping or tearing and the binary 
semi-major axis. Here, we adopt a = 0.1 (r; ~ 50), s = —3/4, q = 0.1 (^ 1^2 = 10/121), 
/3 = '/r/4, and 7 = 7 r/ 4 . The estimated semi-major axes are Owarp ~ 0.01 pc and Otear ~ 
0.02 pc, although they mildly depend on the orbital eccentricity, binary mass ratio, tilt angle, 
and azimuth of tilt, except for (3 = 7 r/ 2 , 0 . Because it is observationally more difficult to 
distinguish between the disk warping and tearing, the actual semi-major axis is likely to be 

®warp ^ ® ^ ®tear’ 

There are observational evidences for disk warping in the maser disks at the center 
of NGC 4258 [10], Circinus [ 6 ], NGC2273, UGC3789, NGG6264, and NGC6323 [9]. The 
tidally-driven warping or tearing of misaligned circumbinary disks can be one of mechanisms 
for explaining the warped structure of the maser disks in these systems. Since no observational 
information is currently available about the warping or tearing radius, we simply assume that 
these maser disks start to be warped or torn at the innermost maser spot radii. Here, we 
pick up NGC 4258 case because of the most remarkable example of disk warping among the 
known maser disk systems. The inner most maser spot radius is 0.17 pc. Assuming that 
-Robs = 0.17pc, the semi-major axis is estimated to be 1.2 X 10“^ pc. However, it is difficult 
to distinguish, solely by the current analysis, whether the central object is a single SMBH or 
such a small-scale SMBH binary. Other independent theoretical and observational approaches 
are needed. 

Probing GWs from individual SMBH binaries with masses > IO^Mq with pulser timing 
arrays (PTAs; [41, 42]) gives a powerful tool to determine if the central object surrounded 
by the warped maser disk is a SMBH binary or a single SMBH. For a typical PTA error box 
(~ 40deg^) in the sky, the number of interloping AGNs are of the order of 10^ for more than 
10® Mq black holes if the redshift range is between 0 and 0.8 (see Figure 1 of Tanaka et al. 
[43] in detail). As discussed in H14a, even if the observed warped maser disk systems host 
SMBH binaries, the characteristic amplitudes of GWs from these SMBH binaries are three to 
four orders of magnitude smaller than the current PTA sensitivity, so that it is unlikely for 
GWs to be detected from these systems. 


4 Conclusion 

We have investigated the tidally driven warping and tearing of a geometrically thin, non-self- 
gravitating circumbinary disk around two SMBHs in a binary on an eccentric orbit, where the 
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Figure 6. Estimate of a semi-major axis of a SMBH binary by the relationship between the 
tidal warp and tearing radii and the observed radins, which is shown by the purple solid line. The 
solid and dashed black lines are, respectively, the tidal-warp and tearing radii of the misaligned 
circumbinary disk around the eccentric SMBH binary with e = 0.6, a = 0.1 (r/ ^ 50), s = —3/4, 
g = 0.1 (C 1 C 2 = 10/121), M = IO^Mq, (3 = 7 r/ 4 , and 7 = 7 r/ 4 . The shaded area between the two blue 
lines represents the whole region of the circumbinary disk. 


original disk plane is misaligned with the binary orbital plane. While tidal torques acting on 
the tilted circumbinary disk tend to align it with the orbital plane, the viscous torque in the 
vertical direction tends to retain its original orientation. The disk can then be warped at the 
tidal warp radius, where the tidal precession timescale equals the vertical viscous timescale. 
On the other hand, the disk is broken into mutually misaligned rings inside the tearing radius, 
where the tidal precession timescale equals the local horizontal viscous timescale. 

There is a critical value of disk viscosity parameter. If a > ^jH/{2>r) with H/r < 0.1, 
the tidal warp radius is larger than the tearing radius. Then, the disk tearing first occurs 
and subsequently disk warping occurs outside the tearing radius. If a < ^^H/{3r) with 
H/r < 0.1, only the disk tearing occurs, because the disk material inside the tearing radius 
is likely to accrete rapidly. 

The tidal-warp and tearing radii most strongly depend on the binary semi-major axis, 
although it mildly depends on the other orbital and disk parameters except for (3 = 0 , tt. 
This strong dependence enables us to estimate the semi-major axis, once the tidal-warp or 
tearing radius is determined observationally. For the tidal-warp or tearing radius of 0.1 pc, for 
instance, SMBH binaries with masses of IO^Mq and other typical orbital and disk parameters 
are estimated to have a binary separation on a 10“^ pc scale. 
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